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ABSTRACT 


It  is  well  known  that  the  Rayleigh  approximation  to  extinction,  scattering  and 
absorption  efficiencies  for  spheres  is  limited  to  small  size  parameters,  z,  and  small  values 
of  |m|x,  where  m  is  the  complex  index  of  refraction.  It  is  also  known  that  the  Thomson 
approximation  to  these  same  efficiencies  is  valid  for  small  x  and  m  -  0 c.  We  have  established 
an  exact  transform  of  the  Mie  coefficients,  for  both  spheres  and  infinite  cylinders,  that 
removes  the  m  related  restrictions  of  the  Rayleigh  and  Thomson  approximations.  The 
resulting  series  for  spheres  and  infinite  cylinders  are  valid  for  all  m  and  small  z.  This 
transform  can  also  be  applied  in  an  approximate  form  to  spheroids. 


RESUME 

L’approximation  de  Rayleigh  pour  1’efficacite  d ’extinction,  de  diffusion  et 
d’absorption  des  spheres  est  liraitee  aux  petites  valeurs  du  parametre  de  taille,  z,  et  aux 
petites  valeurs  de  |m|z,  m  etant  Findice  de  re  fraction  complexe.  L’approximation  de  Thom¬ 
son  est  de  son  cote  valide  pour  les  petites  valeurs  de  z  et  m  =  00.  Nous  avons  etabli  une 
transformation  exacte  des  coefficients  de  Mie  pour  les  spheres  et  les  cylindres  infinis  qui 
fait  disparaitre  ces  restrictions  sur  la  valeur  de  m  dans  les  approximations  de  Thomson 
et  Rayleigh.  La  serie  qui  en  resulte  est  valide  pour  tout  m,  z  demeurant  petit.  Cette 
ti.  .^formation  peut  aussi  s’appliquer  aux  spheroi'des  mais  d’une  maniere  approximative. 
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EXECUTIVE  SUMMARY 

Any  complex  military  system  or  scenario  that  involves  electro-optical  devices 
requires  understanding  the  propagation  of  radiation  through  natural  aerosols  and 
smokes/obscurants.  Given  present  and  future  electro-optical  capabilities  of  fire-control  sys¬ 
tems,  full  spectrum  obscurants  will  be  required  as  effective  countermeasures.  E  is  well 
known  that  spherical  particles  obscure  efficiently  in  the  visible  region  of  the  spectrum  It 
is  also  well  known  that,  in  the  near  and  far  infrared,  as  well  as  the  microwave  region  of 
the  spectrum,  the  use  of  nonspherical  particles  is  mandatory.  The  most  efficient  obscurants 
generally  require  that  the  obscuring  particles  be  made  from  highly  refractive  materials  (far 
example,  metals). 

This  work  considers  the  calculation  of  the  extinction  efficiency  from  small,  non¬ 
spherical  particles  (by  small  particles  we  mean  that  the  particle  size  is  small  relative  to 
the  wavelength).  It  has  two  main  objectives.  The  first  is  to  extend  the  calculation  to 
previously  inaccessible  parameter  domains  and  the  second  is  to  significantly  reduce  the 
computational  burden.  This  would  allow  the  full  exploration  of  the  effects  of  these  particles 
on  the  performance  of  obscurants  and  on  electromagnetic  propagation. 

To  date  the  theoretical  exploration  of  the  effect  of  these  types  of  particles  on  the 
performance  of  obscurants  has  been  either  extremely  restrictive  or  prohibitively  expensive. 
This  work  partially  overt  omes  these  restrictions  and  will  aid  in  the  design  and  performance 
analysis  of  potential  new  obscurants  using  high-refractivity  materials. 

The  electromagnetic  extinction  through  rain,  carbon  smoke  or  other  nonspherical 
atmospheric  aerosols  can  also  be  evaluated  by  this  approach.  This  is  an  essential  first  step 
in  the  evaluation  of  the  performance  of  all  electro- optical  systems. 

The  long-term  goal  of  this  work  is  to  alleviate  the  remaining  constraints  in  the 
theoretical  consideration  of  nonspherical  aerosols  and  obscurants  of  both  natural  and  ar¬ 
tificial  origin.  This  will  not  only  aid  in  finding  better  obscurants  but  also  in  the  possible 
identification  and  remote  classification  of  such  aerosols. 
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1.0  INTRODUCTION 

If  a  particle  is  geometrically  and  optically  small  enough,  then  simple  formulae  can 
usually  be  found  for  the  extinction  and  scattering  efficiencies  as  well  as  the  phase  function. 
Such  approximate  formulae  are  usually  called  the  Rayleigh  approximation  (Refs.  1  and  2).  If 
the  particle  is  still  geometrically  small  but  optically  very  large,  the  Thomson  approximation 
results  (Ref.  3).  Often,  however,  a  geometrically  small  particle  may  be  neither  optically 
small  nor  very  large.  This  creates  a  gap  that  has  not  been  completely  filled. 

In  this  report  we  first  demonstrate  a  transform  that  bridges  this  gap  for  spheres  and 
infinite  cylinders.  If  this  transform  of  the  Mie  coefficients  and  their  series  expansion  were 
limited  to  spheres  and  infinite  cylinders,  it  would  have  only  theoretical  interest.  However, 
this  transform  generalizes,  in  an  approximate  form,  to  other  convex  bodies  when  a  Rayleigh 
or  Thomson-like  series  can  be  obtained.  These  can  sometimes  be  obtained  by  a  method 
given  by  Stevenson  in  1953  (Ref.  4).  Series  are  available  in  this  document  for  ellipsoids 
(which  include  spheroids).  We  will  present  results  of  a  comparison  between  the  newly 
obtained  series  and  the  exact  codes  for  spheres  and  spheroids. 

This  work  was  performed  at  DREV  between  March  and  October  1991  under  PSC 
32A,  EO/IR  Protection  of  Land  Vehicles. 
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2.0  THEORY 


We  will  first  derive  the  general  expansion  of  the  Mie  coefficients  for  small  size 
parameters  i,  and  show  how  the  Rayleigh  and  Thomson  expressions  result. 


2.1  Spheres 


Starting  with  the  definition  (Ref.  2)  of  Q,ca  and  Qext,  the  scattering  and  extinction 
efficiencies,  respectively,  we  have 


Qsca  =  ~2  5Z(2n  +  l){!fln|2  +  l&nl*}, 


Qe*.  =  4  f)(2n+1){R..(r(4  +  &„)}, 


where  an  and  bn  are  the  external  field  Mie  coefficients.  These  Mie  coefficients  are  given  by 


y/fiVniy/V1  X)  Cn(*)  -  y/i  My/W  *)  <»(*)  ’ 


\/f  x)  C»(*)  -  y/&  tM  y/W  x)  Q(x)  ’ 

where  c  and  \i  are  the  relative  dielectric  constant  and  the  relative  magnetic  permeability, 
is  the  Riccati-Bessel  function  of  the  first  kind  and  („  is  the  Riccati-Bessel  function  of 
the  third  kind.  Note  that  the  6„  coefficients  are  symmetrical  with  the  an  coefficients  upon 


substitution  of  c  for  p  and  vice  versa.  Hence,  we  need  only  discuss  an  in  detail. 
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Expanding  aj  and  a 2  in  small  x  we  obtain 


4  (Si) 


3  1  /v-6t  +  4  +  ty\  s 

+  s’  l  (e+2)2  r  + 


iflri  V 

f  +  'll 


1  .  /<V  +  e3(V  +  35p  -  25)  -  e2{70/ii  +  150)  +  200(2f  -  1)\  7 
+  _,  - -  )'  +  ..., 


5n  =  an,  (  *->  fi. 


Expansions  for  higher  order  coefficients  can  also  be  derived.  The  classic  Rayleigh  approxi¬ 
mation  gives  (n  =  1,  m2  =  c/z), 


2  .  /m2  -  1 
Ql  31  1  m2  +  2 


l\  .3.2./  (m2  -  l)(m2  ~  2)\  5  4  /  rn2  -  1 

2/  T5  \  (m2  +  2)2  J  +  9^m2  +  2 


1  .  / (m2  -  l)(m6  -(-  20m4  -  200m2  +  2Q0)\  7 
+  I75I\  (m2  +  2)3  )X  +"-’ 

61  =  7ri(m2  -  l)x5  +  . .  .. 

45 


The  Thompson  approximation  gives  (c  -*  00,  fi  =  0), 


2-3  I.5  4  6  1.7 

ai  —  2*x  +  5**  +  9*  ~Vx 


.  I.3  1.5 

61  =  ~3tX  +  5 12  + 


[10] 
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It  is  immediately  apparent  that  as  m  -*  oo  the  Rayleigh  expansion  for  either  a,  or  6, 
does  not  converge  even  for  small  values  of  x.  Hence,  the  requirement  for  small  optical 
size  jmar|.  The  Thomson  expansion,  having  no  index  dependency,  applies  only  to  infinite 
optical  size.  The  problem  with  the  Rayleigh  series  mathematically  arises  from  the  terms 
in  an  with  t2p  and  in  bn  with  p2e.  If  these  terms  were  set  to  zero,  the  divergence  problem 
would  disappear.  We  must,  however,  recover  the  lost  information  without  reintroducing  the 
divergence  problem. 


To  do  this  rewrite  a„  to  isolate  the  material  properties  as 


an  = 


V’nl>)  -  { 


[11] 


Consider  the  result  of  letting  p  — *  0  in  the  above  expression  (which  will  make  all  the  r2p 
terms  zero  along  with  others), 


lim  an 

p-0 


[12] 


To  recover  the  full  Mie  coefficients,  replace  c  in  the  above  limit  by 


and  by  symmetry, 


\fi  (n  +  1)  tMx/fp  i)  _  s 

^  X  ' 


y/JI(n  +  })xl’n{y/rfix) 

P  -»  "V - 77,  ==— C  =  «n- 


[13] 


[14] 


Both  £n  and  Un  can  be  considered  as  transformed  material  properties.  Putting  these  new 
variables  back  into  [12]  and  expanding  in  small  x,  as  if  £n  andUn  were  independent  of  x. 
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we  obtain, 


2 

Ol  =  X! 


Ui  +2/ 


X3  +  Li(£t-M>+4 


(£1  +  2)2 


Is  + 


f  i  (£>_iiy 

9  Ui  +2/ 


1  /  F,3  +  65?  -  165,  +  8s 


7* 


a2 


(5i  +  2)3 

= 

15*  V252  +  3 ) 


x '  +  .. 


[15] 

[16] 


and  again  by  symmetry, 


bn  =  Qt.,  £«  +-+ ZY„. 


Note  singularities  occur  if  and  only  if 


5n  =  -(n  +  l)/n, 


[17] 


and 


Wn  =  -(n  4-  l)/n. 


[18] 


which  requires  a  real  e  (or  real  //).  For  all  other  values  of  £  and  U  the  coefficients  of  series 
[15]  and  [16]  are  finite  to  all  orders  of  x.  Thus  [15]  and  [16]  are  valid  for  all  nonreal  values  of 
the  refractive  index.  However,  when  the  conditions  [17]  and  [18]  are  nearly  satisfied,  many 
terms  in  [15]  and  [16]  will  be  required  before  the  series  converges  to  a  given  accuracy.  For 
convenience  in  computations  we  can  rewrite  £n  as, 


2  eF  C/A  3c(l+3F) 

£i(z)  =  -77-7  ft(*)=  -7 


(1  +  F) 


2  +  6F-Fz2 


[19] 


with 


F{z)  —  lCOt-(~l — i,  and  z  =  yJTJix  =  (n  -  ik)x. 


[20] 


UNCLASSIFIED 

6 


and 

Since  the  ‘Mie’  coefficients  for  the  normal  incidence  infinite  cylinder  (Ref.  2}  are 
identical  to  the  sphere  Mie  coefficients,  apart  from  the  order  of  t  he  Bessel  functions  involved, 
the  same  procedure  can  be  used  to  obtain  a  series  for  small  particles  independent  of  mate¬ 
rial  properties.  Oblique  incidence  can  also  be  done  but  generates  much  more  complicated 
expressions. 

2.2  Spheroids 


The  above  series  [15]  etc.  by  itself  is  interesting  only  from  the  classic  nature  of 
the  Rayleigh  and  Thomson  series.  However,  the  idea  of  transforming  a  Rayleigh-like  series, 
(with  the  optical  size  constraint)  into  a  series  that  is  independent  of  the  material  properties, 
for  many  types  of  regular  particle  shapes,  would  be  of  practical  use.  We  show  here  that  this 
can  be  achieved,  to  some  approximation,  with  oriented  and  randomly  oriented  spheroids. 


The  Rayleigh  approximation  (with  arbitrary  p)  for  oriented  spheroids  (Refs.  4  and 
5),  Qray ,  is  given  by 

Q  ray  =  Q  sea  +  Qabsi  [21] 


where 


-  !  !£/!!£»  /UU.  x  UM*  x  (I±i2fl£>  (|rf  +  W,»)  },  m 


Q.ca  =  ~  (l^ll2  +  I7?!  I2)  + 


UNCLASSIFIED 

7 


9*  -  Re{.[^  („  + ,/,)  +  <„  + ,/,)  ] }. 


and  where 


p~\Z< cos2  6  +  r2  sin2  8,  a  =  2ttq/A,  6  =  2tri3/A. 


!•)'{? 
I*'*  J 


[24] 


Here  r  =  o/i>  is  the  aspect  ratio  (for  prelates  r  >  1  and  for  oblates  r  <  1),  o  is  the  length  of 
the  semi-axis  of  rotation,  /3  is  the  other  semi-axis  of  the  spheroid,  6  is  the  angle  between  the 
incident  radiation  and  the  a  (or  a  axis),  and  A  is  the  wavelength  of  the  scattered  radiation. 
Furthermore, 


Vi  = 


V2  = 


3(£i  +  ^rr) 


and  r/[  = 


and  t/2 


1 


3(^t  +  7~t) 


3  (£2  +  7^)  ■“  3  (L2  +  ^)* 

and  the  form  factors  are  defined  for  prolates  (i.e.  r  >  1)  as 

(1  ~  92) 


Lx  = 


1  -  Lx 


L*  = 


92  =  1 


r2’ 


For  oblates  (i.e.  r  <  1)  we  have 


Lx  = 


1  +  f2  f  tan  1  / 

P  l  / 

1  -  Lx 


}• 


Lo  = 


2  ’ 


1 


r  = 

r2 


[25] 


[27] 

[28] 
[29] 


[30] 

[33] 


[32] 
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Normally,  p  —  1  in  the  above  equations  and  hence  i/t  and  if2  are  both  zero.  The  procedure, 
as  demonstrated  for  spheres,  requires  the  full  expression  for  the  efficiencies,  therefore  we 
retain  the  p  terms. 

Since  th  j  scattered  wave  at  large  distances  from  finite  convex  bodies  can  be  approx¬ 
imated  by  a  series  of  Riccati- Bessel  functions,  we  will  use  the  same  material  transforms, 
[l3]  and  [14],  but  with  the  arguments  modified.  Furthermore  these  transforms  will  be  angle 
independent  to  first  order.  This  can  be  understood  by  considering  that  the  first  mode  can 
be  excited  in  any  resonant  phenomena  without  regard  to  the  direction  of  the  triggering 
energy.  This  is  re  .dily  verified  by  using  an  >>xact  infinite  cylinder  scattering  code. 

The  modification  is  simply  an  approximation  to  the  effective  optical  size,  z.  of  the 
spheroid.  We  want  a  simple  formula  ’hat  will  give  the  sphere  material  transform  when 
r  =  1  and  gives  the  emipirically  derived  asymptotes  (of  zjb)  as  either  r  —  oc  or  r  — *  0.  It 
can  be  shown  formally  (Refs.  7  and  8}  that  the  rji  and  tf.  terns  come  from  different  vector 
spheroidal  wavefunctions.  This  implies  that  we  will  have,  for  either  prolates  or  oblales.  two 
different  expressions  for  z.  Using  the  T-matrix  method  (Ref.  6)  we  have  empirically  derived 
the  following  expressions: 


y/cju&O  +  v(l  ~  1/r2)) 

prolates 

z  — 

v'cp&(r2r 

oblaies 

>  for  77),  and 

[33] 

v/7/76(l  +  i/1/3(l  -  1  /v'F)) 

prolates 

Z  — 

oblates 

■  for  % 

[34] 

where  v  as  \/3/lG. 
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It  is  found,  when  [33]  and  [34]  are  used  in  the  final  formalism,  the  positions  of  the 


first  two  polaritons  are  well  modelled  from  0.1  <  r  <  20. 


For  randomly  oriented  spheroids  the  cross  sections  C3ca  and  L'a>  -  must  be  integrated 
over  all  angles.  C,ca  is  just  [22]  multiplied  by  irpb2  and  similarly  for  Cab*  from  [23].  The 
efficiencies  are  then  obtained,  in  the  usual  way,  by  normalizing  the  cross  sections  by  the 
average  projected  area.  The  integration  is  simple  since  rp,  tj[,  t}2  and  tjz  are  independent  of 
the  orientation  angle.  Replacing  t  with  £j(z)  and  p  with  U\(z),  integrating  Caca  and  Cab* 
over  sin(0)d0  from  0  to  7r/2  and  then  normalizing  we  get 


~  _  16  64r2  r 

Qsca~  g  ^  {lrl 


(i'll!2  +  \v[\2  2  +  I'Jzl2)}  • 


Qab,  -  3  l*h  +  $1  +  2  [ff2  +  %)]  }  , 


where  the  normalization  factor  is,  for  prolates. 


and  for  oblates, 


A  =  1  +  .  - — sin' 


A=l  +  ~4=^n(Lti^ 
Vl  -  r2  V  r 


^  vr  t  .  j _ ! _ ^  if  r  i  t 


3(Li  + 
1 


3(Li  +  £~r) 

1 


7)2  3(12  +  ^)’  712  3(L2+nd_y  1  J 

Remember  that  in  the  expressions  for  £j  and  U\,  [13]  and  [14],  v/cpr  must  be  replaced  by 


z. 
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3.0  RESULTS  OF  COMPARISON 

In  this  chapter  we  will  demonstrate  the  accuracy  and  utility  of  the  expressions  and 
ideas  in  the  preceding  chapter.  First,  we  will  discuss  the  comparison  of  the  Rayleigh  [7], 
Thomson  [9]  and  our  new  series  [15]  with  the  exact  Mie  solution. 

Figure  1  shows  the  three  series  normalized  by  the  Mie  calculation  as  a  function  of 
size  parameter.  The  series  for  Qaca  are  calculated  to  fourth  order  in  x  and  for  a  refractive 
index  of  500  -  500i  which  corresponds  to  metals  in  the  millimeter  wave  region.  When  [mi] 
is  small  both  the  Rayleigh  and  the  new  series  are  excellent  approximations.  As  expected 
the  Thomson  series  overestimates  Qaca  by  25%.  The  situation  reverses  when  |mij  is  very 
large.  Now  the  Rayleigh  series  underestimates  QtCa  by  25%.  Note,  however,  that  in  the 
intermediate  range  of  jmx|,  [15]  is  still  excellent. 

The  next  figure,  Fig.  2,  shows,  for  the  same  index  as  above,  the  comparison  of  the 
Mie  solution  of  Qext  with  the  new  approximation  and  the  Rayleigh  series.  Now  the  discrep¬ 
ancies  between  Rayleigh  and  Mie  are  quite  obvious.  This  large  discrepancy  is  obviously  due 
to  absorption  (for  more  details,  see  Ref.  1,  p.  292).  This  effect  is  well  modelled  by  the  new 
series. 


To  show  that  [15]  adequately  models  resonances  we  now  choose  an  index  with 
a  large  real  part  and  small  imaginary  or  absorptive  component.  Figure  3  is  a  diagram  of 
Qext  as  calculated  by  Mie  theory  and  the  new  series  [15]  to  fourth  order,  for  an  index  of 
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1.60  -> 


4th  order 


m-500-500i 


1.20  -i 


- Rayleigh 

.  Thomson 

-  New 


I0's  10"*  10'3  1C'*  10''  1 

Size  Parameter  X 

FIGURE  1  -  Three  series  for  Q,ca  normalized  by  the  Mie  calculation,  as  a  function 
of  size  parameter.  Refractive  index  m  =  500  —  500i. 


m=500-500i,  Sphere 


Aporox. 
Rcyleign 
Exact  Mie 


io"  io"  io"  io"  io"  io"  ;o "  io"  i 

Size  Parameter  X 

FIGURE  2  -  Qezt  as  calculated  by  Mie  theory,  the  new  series  [15]  and  Rayleigh. 
Refractive  index  m  =  500  -  500i. 
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Size  Parameter  X 

FIGURE  3  -  Qext  as  calculated  by  Mie  theory  and  the  new  series  [15]  for  an  index 
of  100  -  10_6i. 


FIGURE  4  -  Qtrt  for  water  spheres  at  3.3  GHz.  Refractive  index  m  =  8.743  -  0.6409L 
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Water  9.4  GHz  Sphere 
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FIGURE  5  -  Qext  for  water  spheres  at  9.4  GHz.  Refractive  index  m  =  8.075  -  1 .8241 . 


FIGURE  6  -  Qext  for  water  spheroids,  r  =  2,  at  3.3  GHz.  Refractive  index  m  = 
8.743  -  0.6409i. 
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Water  9.4  GHz.  Aspect  2 
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FIGURE  7  -  (Jejrt  for  water  spheroids,  r  =  2,  at  9.4  GHz.  Refractive  index  m  = 
8.075  -  1.824i. 


100  —  10~6i.  Although  this  index  is  extreme  and  not  physical,  it  is  used  here  to  demonstrate 
that  [15]  is  still  valid  for  such  unusual  cases.  This  emphasizes  that  [15]  applies  to  arbitrary 
indices  for  small  x.  Notice  that  not  all  resonances  are  modelled.  This  is  because  the  higher 
order  Mie  coefficients,  an  and  6n  where  n  >  1,  have  not  been  used  in  this  approximation  of 


Qext- 


The  effect  of  using  the  sixth  order  expansion  instead  of  the  fourth  order  one  can  be 
seen  in  Figs.  4  and  5.  The  cases  are  Qcxi  for  water  at  3.3  GHz  and  9.4  GHz,  respectively. 
It  can  be  seen  that  the  sixth  order  corrects  for  both  amplitude  and  skewness  around  the 
resonances.  Evidently,  higher  order  expansion  would  improve  the  accuracy  further.  It  is 
interesting  that  the  fourth  order  captures  most  of  the  detail,  even  around  the  resonances. 
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Figures  6  and  7  are  the  same  cases  as  for  the  previous  two  diagrams  except  that 
prolate  spheroids,  with  an  aspect  ratio  2,  are  considered  instead  of  spheres.  The  exact 
calculations  are  performed  using  the  T-Matrix  method  (Ref.  6).  Also  the  size  parameter  x 
has  been  replaced  by  the  semi-minor  axis  6.  As  expected,  the  features  are  similar  but  many 
significant  changes  can  be  observed.  Relative  amplitudes  and  locations  of  the  resonances 
have  changed  as  well  as  the  underlying  trend  in  the  curves.  The  agreement,  while  not 
excellent,  is  still  remarkable  because  of  the  simplicity  of  the  series  compared  with  the  exact 
calculation  for  such  a  large  optical  size.  This  large  optical  size  almost  makes  the  T-rnatrix 
ill-conditioned.  No  such  problem  occurs  for  the  series. 

Figures  8  and  9  are  similar  to  Figs.  6  and  7  but  now  oblate  spheroids  of  aspect  0.5 
are  considered.  The  results  are  similar  to  the  prolate  cases  for  both  Qezt  and  the  degree  of 
agreement  between  the  T-Matrix  calculation  and  our  approximation. 

Figure  10  is  a  comparison  between  the  approximation,  T-Matrix  and  the  classical 
Rayleigh.  It  consists  of  a  calculation  of  Qezt  from  a  1.5  aspect  ratio,  randomly  oriented, 
amorphous  carbon,  prolate  spheroid,  at  94  GHz.  This  could  be  potentially  a  part  of  a  mode! 
for  millimeter-wave  propagation  through  a  soot  plume.  Since  both  the  real  and  imaginary 
part  of  the  refractive  index  (m  =  50  -  5Gi)  is  very  high,  the  T-Matrix  calculation  could  be 
unreliable.  Indeed,  for  higher  values  of  6  the  T-Matrix  code  will  produce  negative  values  of 
Qext-  Up  to  about  6  =  0.2  the  code  is  accurate  enough  for  our  purposes  since  the  backscatter 
efficiency  and  absorption  efficiency  behave  as  expected  up  to  this  point  (when  compared 
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Water  3.3  GHz,  Aspect  5 


T-Vctr  * 
Appro* 


*  2  00  4 
<y  5 


t  00  -i 


C  0G 


:  CO  0  2 C  C  *C  0  60  0  60 

Semi-Major  Axis  b 


FIGURE  8  -  Qext  for  water  spheroids,  r  =  0.5,  at  3.3  GHz.  Refractive  index  m 
8.743  -  0.6409i. 


Water  9  4  GHz,  Aspect  .5 

*oc  - 


FIGURE  9  -  Qext  for  water  spheroids,  r  =  0.5,  at  9.4  GHz.  Refractive  index  m 
8.075  -  1.824i. 
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Amorphous  Carbon  94  GHz,  Aspect  1  5 


FIGURE  10  -  Qext  for  amorphous  carbon,  r  =  1.5,  at  94  GHz.  Refractive  index 
m  =  50  —  50i. 


with  equivalent  volume  spheres).  It  is  known  that  these  two  efficiencies  are  more  sensitive 
to  numerical  instabilities  than  the  extinction  efficiency.  Keeping  thi6  in  mind,  the  new 
approximation  agrees  with  the  T-Matrix  calculation  to  a  substantial  degree  while  the  classic 
Rayleigh  series  does  not.  This  is  mainly  due  to  the  inclusion  of  the  magnetic  dipole  terms. 

The  next  example  is  a  model  of  extinction  by  randomly  oriented  copper  flakes  in 
the  infrared  (m  =  35  —  35i).  An  oblate  spheroid  with  an  aspect  ratio  of  0.333  was  used.  For 
this  index,  lower  aspect  ratios  could  not  be  considered  since  the  T-Matrix  will  not  produce 
usable  results  and  hence  no  comparison  could  be  made.  QeTt  for  this  case  is  shown  in  Fig. 
11.  Again,  caution  is  required  since  the  T-Matrix  for  b  >  .25  begins  to  decrease  rapidly 
and  wifl  go  negative,  as  in  the  previous  case.  Again,  as  in  the  previous  figure,  the  new 
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approximation  is  correct  in  both  form  and  magnitude  as  opposed  to  the  Rayleigh  series. 

As  a  fined  example  we  will  use  the  high-temperature  superconductor  YBa2Cu307_£ 
at  77  K.  For  this  material  and  temperature  the  refractive  index  can  be  computed  from  the 
model  presented  in  Ref.  9.  At  10  GHz  we  calculate  rn  =  1.2  —  1 76831  -  The  next  three  figures. 
Figs.  12-14,  show  Qeit  for  a  sphere,  an  aspect  ratio  10  prolate  spheroid  and  an  aspect  ratio 
0.1  oblate  spheroid,  respectively.  As  usual,  the  spheroids  are  randomly  oriented.  In  the 
case  of  the  sphere  the  exact  Mie  calculation  is  given.  For  the  spheroids,  no  exact  solution  is 
known.  The  large  deviations  from  the  Rayleigh  in  all  three  cases  is,  as  in  Fig.  2,  ascribable 
to  absorption.  Notice  that  for  the  largest  x  shown,  Qext  is  close  to  the  value  predicted  by 
the  Rayleigh  series.  Although  it  is  not  clear  from  these  diagrams,  Qcx,  becomes  equal  to 
the  Rayleigh  series,  for  this  range  of  x,  as  the  eccentricity  goes  to  1.  This  is  readily  seen 
from  [39]  and  [40]  since  as  e  — »  1,  L\  — *  0  (prolates)  or  L2  0  (oblates)  and  £\  — ►  oo  and 
U\  — *  0.  Thus  the  electric  dipole  term  dominates  the  magnetic  term,  ^his  dominance  is 
the  weakest  in  the  case  of  spheres  in  which  the  electric  dipole  term  is  only  four  times  the 
value  of  the  magnetic  term  (see  Fig.  1). 

It  is  evident  from  the  above  calculations  that  it  is  possible  to  obtain  series  for  the 
efficiencies  without  consideration  of  optical  size,  for  spheres,  spheroids  and  infinite  cylinders. 
The  accuracy  of  the  series  for  the  spheroid  would  improve  if  the  correct  expressions  for  the 
spheroidal  “Mie”  coefficients  were  used.  However,  the  required  spheroidal  wave  functions 
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Copper,  Infrared,  Aspect  0  333 


Semi-Major  Axis  b 


FIGURE  11  -  Qext  for  copper  oblate  spheroids,  r  =  0.333,  in  the  infrared.  Refractive 
index  m  =  35  —  35i. 


High  Tc  YBCO,  10  GHz,  Sphere 


Size  Parameter  X 


FIGURE  12  -  Qezi  for  spheres  of  YBa2Cu307_i  at  77°K,  r  =  1,  at  10  GHz.  Refrac¬ 
tive  index  m  —  1.2  —  17683i. 
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High  Tc  fBCO,  10  GHz.  Aspect  10 
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FIGURE  13  -  Qext  for  spheroids  of  YBa2Cu3C>7_£  at  77°K,  r  =  10,  at  10  GHz. 
Refractive  index  m  =  1.2  -  176831. 


High  Tc  YBCO,  10  GHz,  Aspect  1 


Semi-Major  Axis  b 

FIGURE  14  -  Qext  for  spheroids  of  YBa2Cu307_s  at  77°K,  r  =  0.1,  at  10  GHz. 
Refractive  index  m  —  1.2  -  176831. 
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would  lead  to  complex  expressions  for  €n  and  Un.  Small-particle  phase  functions  can  also 
be  computed  for  arbitrary  indices  by  using  the  same  material  transformed  series. 

4.0  CONCLUSIONS  AND  REMARKS 

We  have  derived  a  series  that  removes  the  optical  size  constraint  that  is  inherent  in 
both  the  Rayleigh  and  Thomson  approximations.  This  series  has  been  applied  to  spheres 
and,  in  an  approximate  way,  to  spheroids.  Comparison  of  the  new  series  with  the  exact  Mie 
formalism  shows  excellent  agreement.  A  similar  comparison  between  the  T-Matrix  method 
and  the  randomly  oriented  spheroid  series  shows  that  the  agreement  is  not  as  good  as  in 
the  case  for  spheres.  However,  a  substantial  portion  of  the  underlying  structure,  and  hence 
physics,  is  being  modelled  correctly. 

Combining  this  approach  with  our  previously  described  techniques  (Refs.  5  and 
10)  for  approximating  extinction,  scattering  and  absorption  efficiencies,  we  have  obtained 
simple  formulae  for  arbitrary  materials,  particle  sizes  and  aspect  ratios.  For  the  first  time, 
estimations  of  the  extinction  efficiency  for  small,  highly  refractive,  elongated,  spheroidal 
particles  can  be  obtained.  Since  the  phase  function  for  small  particles  is  simply  related 
to  the  first  scattering  terms,  it  can  also  now  be  obtained  for  cases  that  were  previously 
intractable.  This  gives  us  the  capability  to  model  realistic  modern  battlefield  obscurants  as 


well  as  radar  propagation  through  rain. 
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